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We study the topological properties of Peierls transitions in a monovalent Mobius ladder. Along 
the transverse and longitudinal directions of the ladder, there exist plenty Peierls phases corre- 
sponding to various dimerization patterns. Resulted from a special modulation, namely, staggered 
modulation along the longitudinal direction, the ladder system in the insulator phase behaves as a 
"topological insulator", which possesses charged solitons as the gapless edge states existing in the 
gap. Such solitary states promise the dispersionless propagation along the longitudinal direction 
of the ladder system. Intrinsically, these non-trivial edges states originates from the Peierls phases 
boundary, which arises from the non-trivial topological configuration. 
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I. INTRODUCTION 

Many recent efforts have been made to both experi- 
mental and theoretical investigations of the application 
oriented molecular devices^. As new type of the quan- 
tum coherence devices, molecular device emerges various 
novel quantum effects, which enlarge the ranges of the 
material design^i^i^i^. In these systems, the exotic quan- 
tum features would be induced by the non-trivial topol- 
ogy, and the observable quantum effects can also be used 
to specify the topological constructions of the system. In- 
deed, this non-trivial topology induced quantum effects 
never appear in the topologically trivial systems. 

With non-trivial topology, the twisted boundary con- 
dition in the Mobius strip is a good subject to demon- 
strate the significant role of topological structure in low- 
dimensional physics^'''. The Mobius strip is a non- 
orient able manifold, whose edge defines a two-point bun- 
dles over S-*- and thus topological configuration. This 
simple but topologically non-trivial system possesses 
mathematically rigorous description and the accessibili- 
ties of the experimental realization. Actually, the Mobius 
boundary condition has been synthesized in the aromatic 
annulenes, nanographite ribbons and conjugated poly- 
j^gj.g8,9,io,ii,i2^ These progresses motivate us to propose 
a tight-binding quantum device with Mobius topology 
and investigate its quantum properties of the transporta- 
tion of the spinless particles^^'^^. 

Another important phenomenon in molecular de- 



transitions apparently combines both the topological ef- 
fect and the structure instability. In this paper, we 
demonstrate the various dimerization patterns in the 
monovalent Mobius ladder in details. Because of the pos- 
sibilities of the lattice deformations along transverse and 
longitudinal directions, there exist five typical uniform 
dimerization patterns that we will display in this paper. 
In contrast, it is noticed that there is only one dimeriza- 
tion pattern in a one dimensional system. All the five 
dimerization patterns contain the rung, the columnar, 
the staggered dimerization patterns and the vertically 
saw-toothed, the inclined saw-toothed dimerization pat- 
terns as the combinations of the former three ones, all of 
which will be explicitly defined in the next section. 

We compare the Peierls phase diagrams of the Mobius 
ladder with that of the generic one. Here, the generic 
ladder satisfies periodic boundary condition. It is dis- 
covered that when the generic boundary condition was 
replaced by a Mobius one, the conducting properties are 
dramatically changed for the staggered and the inclined 
saw-tooth dimerization patterns and not changed at all 
for the other three patterns. This fact motivates us to 
use the continuum model to analyze the exotic dimeriza- 
tions. We notice the existence of the localized statoi^*^, 
and find the charged solitons propagating in the bulk, 
which promise that the Mobius ladder with staggered 
dimerization is eventually metallic. Being similar to the 
gapless localized states in graphene strip, we also point 
out that our model behaves like a "topological insula- 

,7721.22 



vice^ii^ is their Peierls instability, which exists univer- tor"i^i^ with localized state existing at topological non- 



sally in low-dimensional physical system including poly- 
mers, spin chains, and organic materials, etoi^^. The 
significance of the Peierls transition is that after the lat- 
tice is deformed due to the electron-phonon coupling, the 
system energy is decreased and thus the changed energy 
band structure converts the original metal phase into an 
insulator one. Actually, the existence of metal-insulator 
transition in the polyacetylen o^^i^^ originates from this 
lattice modulation. 

For the Mobius ladder configuration as a quasi-one 
dimensional (QID) system, investigation of the Peierls 



trivial boundary. Here, the topological insulator refers to 
a bulk insulator which possesses robust metallic localized 
states, which is different from mundane band insulator. 
These localized states are actually topologically in- 
variant, which characterize the time-reversal invariance 
of the topological insulator. 

This paper is organized as following. In Sec. [Ill we 
present the lattice Hamiltonian of the Mobius ladder and 
calculate the Peierls phase diagram. In comparison with 
the Mobius case, we also calculate the Peierls phase di- 
agram of the generic ladder. To prove the existence of 
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FIG. 1: (Color online) Schematic illustrations of the ladder 
with (a) Mobius boundary condition, (b) generic boundary 
condition and (c) the corresponding one-dimensional version 
of Mobius ladder system with long range coupling. 

the non-trivial localized states, we introduce the contin- 
uum model of the Mobius ladder in Sec.[TTllas well as the 
one of the generic ladder without any solitonary solution. 
We conclude our main results in Sec. [iVl The detailed 
derivation of the continuum model from the lattice model 
is shown in Appendix „ 

II. TOPOLOGICAL PEIERLS TRANSITIONS 
AND CORRESPONDING PHASE DIAGRAM 

A. Model setup and dimerization patterns 

In this section, we describe a tight binding model for 
the electrons hopping on a ladder with the tight-binding 
Hamiltonian 

N-l N-1 

= E A]m, A, - E J, A] A,+i + /..c, (1) 

where operator- value vectors Aj = (aj^bj)^ is defined 
in terms of the annihilation operator aj and bj of the 
upper and the lower chain of the ladder (see Fig. [TJa)), 
which are denoted as a-chain and 6-chain in the following 
discussion. The transition matrices 

= sja, - Gjaa: (2) 

and Jj = diag[Jj,Jj^] are defined by Pauli matrices 
(Jx^cTy and cr^, on-site energy differences Sj = £o, cou- 
pling strength between a-chain and 6-chain Gj = Go and 
hopping strength Jj = = Jq. Here, N is the site num- 
ber of the a- (6-) chain. To demonstrate the effect of the 
topology configuration of such system, we will consider 
two kinds of boundary conditions, which are the Mobius 
boundary condition characterized as 

Aj+AT = CFx-^j (3) 
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FIG. 2: (Color online) Schematic illustrations of electron fill- 
ing for the monovalent ladder system before Peierls transi- 
tions, which include (a) insulating case and (b) metallic case. 

(see Fig. [Tfa)) and the generic one characterized as 

Aj+AT = A^- (4) 

(see Fig.mb)). 

The boundary condition apparently affects the sys- 
tem globally, and different boundary conditions result 
in different symmetries. The generic boundary condi- 
tion corresponds to rotational symmetry, which implies 
the generic ladder possesses topological configuration. 
In contrast, the Mobius ladder is considered as a non- 
orient able manifold, whose edge defines a two-point bun- 
dle over and thus topological configuration^^. This 
unusual topology can induce some novel effects such as 
induced gauge field and the cut-off of the electrons trans- 
mission spectrum^^. In our paper, the topological 
configuration will contribute to the formation of the non- 
trivial localized states. 

By diagonalizing the above tight binding model ([T]), 
two energy bands 

Ek = ±^el^Gl + Jo cos (k + ^) (5) 

can be obtained for the Mobius ladder, and for the generic 
ladder the corresponding two energy bands are 

Ek = ± ^jel^Gl + Jo cos k. (6) 

The quantity of the energy shift due to the ii/N phase 
shift in the energy spectrum (Eq. ([5|)) of the Mobius 
ladder depends on the position of the level and results 
from its nontrivial topologyi^. If we only consider the 
monovalent case that 2N electrons are filled in the all 
negative levels for the ladder system, there are only two 
kinds of energy spectra and the corresponding filling con- 
figurations (see Fig. [2]) before the Peierls transitions. 
One case is that the valence band is entirely filled by 
the elect rons and the conduction band is empty when 
2Jo < v^o+^O' which corresponds to the insulator or 
the semi-conductor phase (Fig. EJa)). Another case is 
that the electrons fill part of the conduction band when 
2Jo > yio+^O' which corresponds to the metal phase 
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(Fig. [2](b)). Since the Peierls transitions discussed below 
actually change the phases of the system from conductor 
to insulator, only the second case is taken into account 
in the following discussion. 

In order to consider Peierls transition induced by 
electron-phonon interaction in the ladder system, we use 
the Born-Oppenheimer approximation by presuming the 
transverse and the longitudinal lattice deformation (see 
Fig. [3]) depicted by two displacements 6 and a, which 
are small comparing with the lattice constant of the 
transverse direction / and the longitudinal direction m. 
In Fig. [3l we have assumed that the lattice is uniquely 
dimerized, where in general cases the deformations de- 
pend on the locations of the sites. The above approxi- 
mation is valid because the frequency of the phonon is 
much smaller than the frequency of the electrons. In the 
sense of the Born-Oppenheimer approximation, we can 
fix the displacements of the lattice to solve the eigenval- 
ues of the electrons, which eventually act as the effective 
potential onto the phonons. 

Since the transition matrices Mj, Jj of the electrons 
in Eq. ([T]) depends on the distance between the nearest 
neighbour sites, after the lattice deformation the transi- 
tion matrices depend on the displacements S and a as 
well. Additionally, the lattice deformation is modeled as 
2N coupled harmonic oscillators with the Hamiltonian 
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FIG. 3: (Color online) Schematic illustration of the dimeriza- 
tion patterns including (a) original lattice, (b) transverse, (c) 
columnar, (d) staggered, (e) vertical saw-tooth, and (f) in- 
clined saw-tooth, m and / are the lengths of transverse and 
longitudinal directions, respectively. 6 and a denote the static 
deformations along transverse and longitudinal directions, re- 
spectively. 



clined saw-tooth (Fig.[3Kf)), are regarded as that the sys- 
tem is transversely and longitudinally dimerized simul- 
taneously, which possess either columnar dimerization or 
staggered dimerization along longitudinal direction. 



^(ma,j - rubj)'^ + X] X] ^ (^^ j) X^) B. Band spectral structures of the Mobius ladder 



J=0 



i=a,b j=0 



where kj and rriij {i = a, b) are the displacements of the 
j-th site of a-chain or 6-chain along the longitudinal and 
the transverse direction, respectively. Here, Kt and Ki 
are spring constants of the transverse and longitudinal 
directions, respectively. 

The Peierls transition happens when the decrement of 
the total electrons energy compensates the increment of 
phonon energy caused by the lattice deformation. Here, 
the Fermi surface plays an important role. Actually, after 
the lattice is modulated, the gaps, which are opened up 
at the Fermi surface, mainly result in the decrement of 
the total electrons energy. 

All five uniform dimerization patterns , including three 
simple dimerization patterns and two hybrid ones, are 
presumed for manovalent case, which is illustrated by a 
deformed two-square section of the ladder in Fig. [31 The 
undimerized lattice is present in Fig. [31(a). The first sim- 
ple case is the rung dimerization (Fig. [S^b)) which pos- 
sesses lattice deformation only along the transverse direc- 
tion. Along the longitudinal direction there are two dif- 
ferent dimerization patterns: the columnar dimerization 
(Fig. ^c)) and the staggered dimerization (Fig. [S^d)), 
which correspond to same or different Peierls phases in 
the a— and the b— chain. Last two hybrid dimerization 
patterns, the vertical saw-tooth (Fig. [SJe)) and the in- 



Since the Fourier transformation is no longer valid 
when the boundary is not periodic, we introduce a 
new operator- value vectors Bj = UjAj, where the site- 
dependent unitary transformation^^ is defined as 



exj[){—iOjl) — ex.p{—iOjl) 
1 1 



(8) 



where 6 = tt/ {Nl) is half of the momentum quanta. 
Through this transformation, the period boundary con- 
dition Bj+AT = Bj is retrieved for the new operator- value 
vectors. In the new representation, the Hamiltonian 

N-1 N-1 

He=Y. b1m;.B, - j;BtB,+i + h.c. (9) 

is unitarily transformed from Eq. ([1]), where the new 
transition matrices 



Sj ex.p{iOjl) 



Sj ex.p{—iOjl) 



and 



J- 







(10) 



(11) 
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differ from the original ones. Such difference is consid- 
ered as the global effect induced by the Mobius boundary 
condition, where an induced gauge field takes responsi- 
bility for the cutoff of the transmission spectrum and a 
stark shift occurs in the energy spectrum^^. However, 
this two effects actually are not significant when we only 
consider electrons filling in the energy bands for static 
dimerization with very large site number. 

To obtain the static lattice deformations, it is neces- 
sary to minimize the total energy E = E^^E^ versus the 
lattice deformation. Here, based on Born-Oppenheimer 
approximation, the total energy including the electron 
part Eq and the phonon part is obtained by diagonal- 
izing the electron Hamiltonian Eq. (|9]) and the phonon 
Hamiltonian Eq. ([7j), respectively. 

We take the staggered dimeriation of the Mobius lad- 
der as an example. Let the m and / be the lattice con- 
stants along the transverse and the longitudinal direc- 
tions, respectively, and we can define the static uniform 
deformations Uj = { — ly 6 and Vj = { — ly a. Because 
the lattice deformation changes the coupling strength and 
the hopping strength from Go and Jo to Go-\- (3{m^ — m^) 
and Jo + a{lj — /j+i), where a and /3 are the rate of the 
changes of the longitudinal and the transverse hopping. 
For the staggered dimeriation, the longitudinal lattice de- 
formation is Ij = I -\- {—iy~^^ a and the transverse lattice 
deformation is 

- = \/^^ + (2cr)^ - m ^ 2crVm. (12) 

Therefore, with the modified coupling strength Gj = 
Gq = Go -\- 2Pa'^/m and hopping strength Jj = Jo + 

AJ(-l)^' for a-chain and = Jo + AJ(-1)^'+^ for b- 
chain with A J = 2a<j, four separate energy bands of the 
electrons can be obtained by diagonalizing the Hamilto- 
nian Eq. (|9]) in the momentum space as 

ej{k) = {-iy^J^{a) + {-l)li\i.ia), (13a) 
^(ct) = Gl+ A J2 + 4J2 cos2 (kl + ^) , (13b) 

Kct) = 2^G'o^AJ^ + 4JiGl cos2 (^kl + ^) (13c) 

for j = 1,2,3,4, where , and [|J represents the integer 
part of of j/2. 

It follows the energy band diagram (Fig. S^b)) that 
when the hopping strength Jo is sufficiently large, the 
deformation opens four gaps in the original two over- 
lapped bands. The two gaps at /c = ±7r/2/ are usual 
ones because they only arise from the longitudinal defor- 
mation for a- and 6-chain, respectively. The other two 
gaps approximately locating at Fermi momentum kf = 
± arccos(Gm/2Jo)/^ in the upper band and k'^ = n /l — kf 
in the lower band basically arise from the coupling be- 
tween the k— states in a— chain and the /c — tt//— states in 
6-chain with strength aaG'^j^Jo approximately. It is es- 
sential to indicate that when the the hopping strength 
Jo is sufficiently small, there are no gaps opened up 
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FIG. 4: (Color online) Schematic spectra of ladders (a) with- 
out dimerization (a) and (b),(c) with staggered dimerization. 
In (b), four gaps are opened up at Fermi momentum /c/, k'^ 
and k — while in (c) two gaps up at Fermi momentum 

disappear when the hopping strength Jo is sufficiently small. 
The shadow regions represent the electron occupation in the 
energy bands. 

at Fermi surface even there exists the coupling between 
the /c-states in a-chain and the k — 7r//-states in 6-chain 
(Fig. SJc)). Thus no Peierls transitions occur. 

To compare with the result of the Mobius ladder, we 
also consider the dimerization in a generic ladder shown 
in Fig. [T](b). Here, the Fourier transformation is applied 
to diagonalize the electron Hamiltonian without intro- 
ducing site-dependent unitary transformation. The stag- 
gered dimerization is still taken as an example for the 
generic ladder, which also possesses four separated en- 
ergy bands 

e'j{k) = (-1)Vm'(^) + (-1)L'J^'(^), (14a) 
/i'(cr) = + A J2 + 4 Jo^ cos^ (^0 . (14b) 
iy\a) = 2^£2^J2+4J2G^ cos^ {kl), (14c) 

which also follows the energy band diagram Fig. S^b) 
and Fig. [D^c). However, the gaps opened up at Fermi 
surface are different from the ones of Mobius case, which 
results in different energy of each dimerization patterns 
of generic ladder from the ones of Mobius case and thus 
the different phase diagrams. 

C. Phase diagram of the Mobius ladder and the 
generic ladder 

Now we focus on the case that the gaps opened up at 
Fermi surface may decrease the energy of the electrons 
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by AEe = E{5) - E{0), where 



= 1,3 



(15) 



and £^(0) is the energy without dimerization. The lattice 
deformation also increases the energy of phonons by 



AE^ = AKiNa'^ + 2KtN- 



m 



2 • 



(16) 



The total energy shift AE = AE^ + AEp versus lattice 
deformation a is plotted in Fig. [5l There is a minimum 
of AE at (T = cTg {Ki^ Kt), which corresponds to the sta- 
ble configuration of the system. Obviously, as the or- 
der parameter of the staggered Peierls phase transition 
in Mobius ladder, as (Ki^Kt) = means the lattice is 
not deformed corresponding to the original metal phase, 
while as (Ki^Kt) ^ means the lattice is spontaneously 
modulated to form an insulator phase. In this sense, 
when all the possible Kt and Ki are chosen to determine 
respective stable configurations, we obtain the phase dia- 
gram of the staggered dimerization in the Mobius ladder. 

The above calculation is carried out for the staggered 
case. Repeating it for all deformations (Fig. [2]) gives 
the total Peierls phase diagram for the Mobius boundary 
condition (Fig. [SJa)). Here, the parameters are chosen 
as Go = 15so, -^0 = lO^o, a = P = eo/m, and / = m. 

With these parameters, only three dimerization pat- 
terns survive, which are the rung, the staggered, and the 
inclined saw-toothed dimerization patterns (all are de- 
noted by capital letters "R", "S", and "I" in Fig.[5Ka)). 
The three Peierls phases occur at different regions at Ki 
and Kf. The rung dimerization occurs when Ki ^ Kt^ 
the staggered dimerization occurs when Ki <C Kt^ and 
the inclined saw-toothed dimerization occurs when Ki ^ 
Kt is sufficiently small. When Ki ^ Kt is sufficiently 
large, there is no dimerization emerging in the Mobius 
ladder. 

With the same procedure, by minimizing AE of the 
generic ladder, we obtain the total Peierls phase diagram 
for the generic boundary condition (Fig. [5](b)). Notice 
here that the energy bands completely filled with elec- 
trons become e'j{k) {j = 1,3) in Eq. (|14ap . 

With the same parameters of the Mobius case, the ba- 
sic properties of the Peierls phases of generic ladder is 
similar to the one of Mobius case. However, the region 
of the staggered dimerization pattern under the Mobius 
boundary condition shrinks comparing with the generic 
one. This fact means that the metal phase is prefer- 
able for a Mobius ladder system. Therefore, the above 
phase diagrams show that the conducting properties can 
be dramatically changed in when the topology of the lad- 
der is switched. Because the inclined saw-toothed phase 
contains the staggered dimerization along longitudinal 
direction, it is changed the same way as the staggered 
one. The Peierls phase for rung dimerization is exactly 
the same whatever the boundary condition is. Although 
the columnar and the vertical saw-toothed dimerization 
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FIG. 5: (Color online) The phase diagrams of the ladder sys- 
tem with (a) the generic and (b) the Mobius boundary condi- 
tions, which are plotted versus {Ki,Kt) in (a) and (b), respec- 
tively. The parameters are chosen as Go = 15so, -^o = lOso, 
a = P = eo/m, and / = m.The distribution of the total en- 
ergy versus {Ki,Kt) determines the boundaries of the 
phases, which are plotted as dashed lines. Here S, I, R, and 
N represent the staggered, the inclined saw-tooth, rung and 
no dimerization, respectively. 



do not appear in the current phase diagrams in Fig. O 
their corresponding Peierls phases are the same for dif- 
ferent boundary conditions of the ladder system when 
Go <C Sq. Therefore, to further consider the topologi- 
cal effect on conducting properties, only the staggered 
dimerization pattern is taken into account and the exis- 
tence of the localized states will be revisited for our QID 
system. 



III. CONTINUUM MODEL AND SOLITONARY 
SOLUTIONS 

A. localized states in the Mobius ladder 

We adapt the continuous field approach by regarding 
the Mobius ladder as a one-dimensional system with long 
range hopping (Fig.[Hc)). The detailed derivation of the 
continuum model from the lattice Hamiltonian Eq. ([1]) 
and (0 is presented in App. „ Without loss of the gen- 
erality, we focus on the special case with eq = 0. For a 
continuous field approach, it is crucial to introduce an 
order parameter 



A{x) = -4a(/)(x), 



(17) 



where (j){x) is the continuous limit of = {—lY Uj. 

The Hamiltonian of the continuum model H = He+Hp 
contains the the phonon part 



M 



[A(a; + L)-A(x)]"| 



(18) 
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and the electron part 



He = J dx<^>{x)Jf^<^>{x), 



(19) 



where M is the mass of the particle and L = Nl is the 
length of the a-chain (6-chain), which approaches infin- 
ity at the end of the calculation. Here, the subindex m 
stands for the Mobius ladder. In the electron part, the 
hopping electron could be described with a 4-component 
spinor <l>(x) = [ (j)2{x) 03(^) 04(^) ] • Physi- 

cally, (j)i{x){(j)s{x)) and ^2(^)(^4(^)) respectively repre- 
sent the left-traveling wave and right-traveling wave in a- 
chain (6-chain). In this spinor representation, the Hamil- 
tonian density is expressed by Pauli matrices a^;, a^, 
as 



^ — 



mazdx + A{x)ax G (x) 

G (x) mcjzdx + A(x + L)cFx 



where = 2/ Jq, crz-,crx are Pauli matrices, and 



G(x) =Go- 



^ [A(x)-A(x + L)]' 



(20) 



(21) 



is effective coupling between the a-chain and 6-chain. 

To refiect the Mobius boundary condition in our QID 
model (Fig. [U^c)), we take the period 2L for boundary 
conditions 



^{x^2L) = ^ {x) 



(22) 



rather than L for the generic case. With this bound- 
ary condition, we solve the Bogoliubov-de Gennes (BdG) 
equation 



^m^i {X) = Si^i {X) , 



(23) 



where i represents the i-th energy band of the spectrum 
and {x) = [(pl {x) , 0| {x) , 0| (x) , 0i (x)]^. At zero 
temperature, the order parameter A {x) satisfies the self- 
consistent equations 



-Afx) - 



{A{x ^ L) - A{x)f 




(^) ^2 (^) 
</>3* (^) ^4 (^) 



, for X < 0, 
, for X > 0, 



(24) 



which is obtained by the functional variation of 

E{A{x))=Y,£i + Hp (25) 

i 

with respect to 6 A (x) and 6 A {x -\- L) . The sum is over 
the energy levels below the Fermi surface. In principle, 
the eigenenergies e^, the eigenfunctions (x) and the 
order parameters A(x) and A (x + L) can be completely 
determined by the BdG equation in Eq. ([23|) and the 
self-consistent equation in in Eq. ([23 



After introducing new 4-component spinor ^i(x) = 

[(p{{x) ifiix) ifiix) (fiix)]^ by ^i{x) = U^i{x), 
where the unitary matrix is defined as 



U 



1 

V2 



1 i 
1 -i 
1 z 
1 -i 



(26) 



the BdG equation in Eq. (|23|) can be simplified as 

^^^,(x) =s,$,(x) (27) 



with new Hamiltonian density 



mcjxdx - A{x)(jy 
Go 



ma^dx 



Go 

- A(x ■ 



L)(Ty 



(28) 

As we shown as follows, some solutions of the above 
BdG equation can exist as localized states. In the fol- 
lowing, we only consider the case Ki <C i^^. In this case, 
three dimerization patterns of rung (Fig. [Sjb)), vertical- 
saw tooth (Fig. [Sje)) and inclined saw-tooth (Fig. [S^f)) 
occur rarely. Thus we only need to compare the energy 
of staggered dimerization with columnar one. Here, we 
revisit Mobius ladder system with the staggered dimer- 
ization characterized by 



A (x) = - A (x + L) . 



(29) 



In this phase, the order parameters in a-chain and b- 
chains are opposite and display a Peierls phases domain 
wall when the site number N is even. We also notice that 
the columnar dimerization to be compared is character- 
ized by 



A (x) = A (x + L) . 



(30) 



For the staggered case, we assume a kink deformation as 
the form 



A(x) = Atanh(x/^) 



(31) 



with ^ = Vf/A, which is so small that the effective cou- 
pling between the a-chain and 6-chain G (x) ^ Go- 

To solve the new BdG equation, some symmetries of 
the Hamiltonian density can be used to simplify the cal- 
culation. The Hamiltonian density actually possesses the 
discrete symmetry 



with an anti-diagonal matrix 



1 
10 
10 
10 



(32) 



(33) 



denoting the mirror refiection transformation. This sym- 
metry guarantees that if ^i(x) is an eigen function of 
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with eigenenergy Si, the W^'^i{x) is also the eigen func- 
tion of with the same eigenenergy Si. Obviously, the 
eigenvalues of matrix W are ±1. Together with the trans- 
lational symmetry characterized by momentum quantum 
numbers, the total Hilbert space can be spanned by these 
bases. Therefore, the eigen function ^i(x) either has the 
form 



(34) 



-kt-kc 



(a) 



kskt 



C 



(b) 



or 



(35) 



In this sense, as the solutions of the BdG equation 
with energy £s = 0, two degenerate solitonary states can 
be found as one with non- vanishing components 



^\ {x)=ipl {x)=FX{x), 
and another with non- vanishing components 
^l{x) = -^l {x)=Ft{x) 

for 



Fi{x) 



\ exp I ±z— X I sech , ^ 



(36) 



(37) 



(38) 



where the subindex s denotes solitonary solutions. These 
solitonary states are the localized states located at the 
midgap. Since there is no such solitonary state in the 
generic ladder, the existence of the solitons is absolutely 
topological effect. 

We note that the another two bands 



{vik ± GoY + A 



(39) 



(illustrated in Fig. (Ua) as two overlapped shadowed do- 
mains) fully occupied by the electrons correspond to 
eigen functions 



^1 (x) 



±(^I'±(x) = 



2/L 
1 

2/L 



-ikx 



Fl{x)e- 



-ikx 



where 



Flix) 



■ Ax .{vik±Go) 

tann h i- 

Vf A 



(40a) 
(40b) 

(41) 



represents a deviation from a plane wave in the kink or- 
der. Here, the subindex v stands for the valence bands 
below the Fermi surface. 

Then, it follows from the self-consistent equations Eq. 
that 



2Kf 



Ki 



tanh^(x/0 



l_ n 



dk 



(42) 



FIG. 6: (color online) Schematics of the energy spectra 
of the valence bands under (a) staggered and (b) colum- 
nar dimerizations, where — =F Gq/v^ and k^ — 



Gl + vjkl ± 2Go^MWTA^/vi. The shadow regions rep- 
resent the electron occupation in the energy bands, and the 
brown straight line represents the solitonary states. 



Usually, the solitons are localized around the original 
point with the width ^ of several lattice constants, which 
is much smaller than the total length of the ladder sys- 
tem. For most sites far away from the original point as 
kl ^ 1^1 5 we assume tanh^(x/^) ^ 1. Since the second 
term is small comparing with the first term at the left 
side of the above equation, we obtain approximate order 
parameter 



Ao exp(-5Ao 



(43) 



where Aq = W exjp{—A) is the order parameter for one 
dimensional uniformly dimerized system, W = 2vfkf, 
A = VfirKi/Sa'^l, and B = Vf7rKt/2^m'^a^l with /cf is 
the Fermi momentum. 

Moreover, we can further prove that the above soli- 
tonary states are the ground states. To this end we 
calculate the total energy of the electron-phonon system 
according to the phase shift^^ of the eigenstates. The 
phase shifts are determined by the eigenstates of the band 
electrons when x tends to ±oo as 



lim F^{x) — -rr 



(^f ^ ± Go) 



(X exp(zi9|^ (/c)). 



which reads 



arctan 



v^k ± Go 



'vfk±Go\ 
arctan : — tt. 



(44) 

(45a) 
(45b) 



Therefore, the total phase shift of the eigenstates is de- 
fined by their difference as 



0^ {k) = 0^^ {k) - etoo (A:) = TT + 2 arctan 



Vfk ± Go' 



(46) 
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A straightforward algebra explicitly give; 



,18 



TT Vfkf 



(47) 



where E^ is the total energy for the columnar dimeriza- 
tion (Fig. He)), and S{Kt) = 13KtA^Vf/{3 x A'^m^aH) 
results from the coupling between the a- and 6-chain. 
The second term in E^ is usual energy increment due 
to the existence of solitonary states. The third term 
in E^ results from the difference in the total ener- 
gies in two filling ways. One corresponds to the stag- 
gered dimerization (Fig. [D^a)) with two lower bands 



— Y (vf/c ± Go)^ + occupied by electrons, while the 
other corresponds to columnar one (Fig. (Ub)) with two 



lower bands ±Go — y {v^k)'^ + occupied. For the lat- 
ter the energy of electrons increases because a part of 
electrons are forced to occupy higher energy levels. If Go 
is so large that dE — E^ — E^ is negative, the ground 
state of the Mobius ladder system corresponds to the 
staggered dimerization rather than the columnar one. In 
this case the solitonary states are localized states as the 
ground state. 



B. Comparison with the generic ladder 

It is a complete topological effect that the ground state 
is localized. To demonstrate this, the continuum model 
for the generic ladder H' = -\- is presented to 
compare with the Mobius case. Here, the phonon part is 



E 

c=a,b 

f 

J — ( 



dx- 



[Aa(x) - Ab(x) 



and the electron part is 

POO 



(48) 



(49) 



{x + L) = (x) with period L. We will show that be- 
cause of the trivial topology of the generic ladder, there is 
no solitonary solution for the dimerization of the generic 
ladder. 

At zero temperature, the order parameter A (x) satis- 
fies the self-consistent equations 



4aH ' A^w?aH 
^2Re (x))*4''(x) 



{K{x) - Ab(x))3 



4aH 
^2Re 



(Aa(x) - Ab(x))3 



(52a) 



(52b) 



respectively, which are obtained by the functional varia- 
tion of 



i;(A(x)) = ^£i + ffp 



(53) 



with respect to (5Aa {x) and (5Ab {x) . The sum is over the 
energy levels below the Fermi surface. For the staggered 
case, we assume a kink deformation Aa {x) = — Ab {x) = 
A. 

After applying the same transformation, the BdG 
equation of the staggered dimerized generic ladder is 



JP^^ix) = e\%{x) 



with new Hamiltonian density 



mcTxdx - Aay 



(54) 



(55) 



where G[) = Go + pA^/^ma^. 

The exact spectrum solved from the BdG equations 
contains four energy bands including the lower two bands 

e'^^ = —\J {vfk ± Go)^ + A occupied by the electrons 

and the higher two bands s'^^ = \J {vfk ± Go)^ + A with- 
out being occupied. The corresponding eigen function are 



where the 4-component spinor (x) has the same physi- 
cal meaning as the one in the Mobius case, and the Hamil- 
tonian density reads 



^ = 

where 



ivfa^dx + Aa(a;)(Tj; 



G{x)=Go 



G'q 

muzOx + Ab(a;)cra: 



[Aa(a;) - Ab(a;) 



(50) 



(51) 



is effective coupling between the a-chain and 6-chain. 
It is noticed that the order parameters here are no 
longer unified in one chain, and are defined as Aa(x) 
and Ab {x) for a-chain and 6-chain, respectively. Addi- 
tionally, the boundary condition for the generic ladder is 



v,± 



1 i(^-k±et 



2/L 



(56a) 



^I'g^ (x) = (x) = , (56b) 



2/L 



for the two lower bands and 



(57a) 



(^i;^ {x) = ±(^4;^ {x)=±^e V 

^2,g (^) = ±^3,8 (^) = ^y^e 



(57b) 



for the two upper bands, respectively. Here, the subindex 
g stands for the generic ladder. 
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There is no solitonary solution existing for the stag- 
gered dimerization pattern of the generic ladder. There- 
fore, the localized states are the complete topological ef- 
fect. 

The Mobius ladder with staggered dimerization ac- 
tually behaves like a topological insulator. Naturally, 
the Mobius configuration is 1? topologically invariant 
and gapless localized states exist in the gap. Adition- 
ally, the topology of the system can protect the soli- 
tonary states from external perturbations. For example, 
when the soliton propagates along the longitudinal di- 
rections without spreading, the energy increment caused 
by moving soliton with velocity from the time evolu- 
tion of order parameter A(x,t) = Atanh(x — v^t)l^ is 
A^^s = A^/(24vfa^/), which could be much smaller 
than the exciting energy bE. It indicates that the mov- 
ing solitons can propagate in the Mobius ladder without 
dispersion and thus is robust to external perturbations. 



IV. CONCLUSION 

We study the topological properties of Peierls transi- 
tions in a monovalent Mobius ladder in contrast to the 
Peierls transitions in a generic ladder. According to lat- 
tice deformation along the transverse and longitudinal 
directions of the ladder configuration, there exist plenty 
Peierls phases corresponding to various dimerization pat- 
terns. The insulator phase resulted from staggered mod- 
ulation along longitudinal direction behaves as a topo- 
logical insulator, which is different from mundane band 
insulator. Actually, this non-trivial insulator originates 
from the Peierls phases boundary induced by the non- 
trivial 1? topological configuration. 
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APPENDIX: THE DERIVATION OF THE 
CONTINUUM MODEL 

When the site number N is so large that the charac- 
teristic wave length of the eigen function is greater than 
the lattice constant /, the continuum field approach is 
appropriately adapted by regarding the Mobius ladder 
as a one-dimensional system with long range hopping 
(Figlljc)). We consider the upper chain (a-chain) and 
lower chain (6-chain) as the first half and second half of 
a whole chain with 2N sites, which corresponds to the 
mapping 



A 



■j+N 



(A.l) 



with fermionic operators Aj. After the mapping, the elec- 
tron part of the lattice Hamiltonian in Eq. ([T]) with a 
fixed deformation configuration {uj} = {uj^u^j} can be 
rewritten as 

N-l N-1 



2N-1 



j=0 



2N-1 



j=N j=N 



where the modified coupling constants are 
G 



^ " ' 2m 

Jj = Jo ^ a (i^j+i - u'^j) , (c = a, b) . 



(A.2) 



(A.3a) 
(A.3b) 



and indecies a, b stand for the original a-chain and b- 
chain. The energy of the phonon in Eq. Q can also be 
obtained as 



N-l 



N-l 



J=0 
2Ar-i 



Ef (»?) 

2 '^'M 



j=N 



j=N 



E 

3=0 



Am? 



(A.4) 



Usually the wavefunction varies greatly from site to 
site under dimerization, which means the coordinate j is 
not suitable for the continuous field approach. However, 
if we introduce the new coordinate 



(2j + ^K 



(A.5) 



as the center of the 2j-th and (2j + l)-th sites, the wave- 
function varies slowly and the continuous field approach 
is valid. 

In this sense, the new fermionic field operators 



(A.6) 



{^l {Xj) = i-iy {iA2j + A2j + i) , 
(Xj) = ^ (-1)^' {A2J + iA2j+i) , 

which satisfy the anti-commutate relations 

|(/Pc(^j)^^d(^jO}^ = ^c,d^^j,/,(c,d = l,2XA.7a) 
{(Pc{xj) ,(pd{xj')}^ = 0, (A. 7b) 

corresponds to the slowly varying wavefunctions of the 
new coordinate Xj. Thus the field operators at Xj can be 
expanded as 



(fc (Xj^l) = (fc (Xj) + 



d(pc (x) 



dx 



21 



,(c = l,2). 
(A.8) 
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The dimerization implies the deformations = 
(— 1)-^ (/)^, which leads the displacement order parameters 

= -2a((/)^2,+</>2,+i)- (^-9) 

As 2N is very large, the summation can be substituted 
with the integral 




(A.IO) 



as well as the field operators 

ipc{xj) ^ V2lipc{x). (A.ll) 

Obviously, from Eq. ()A.7aP and ()A.7bp . the above field 
operators cpc (x) satisfy the anti-commutate relations 

{(^c(^),^i(^0}^ = (^c,d^(x-xO,(c,d = 0(,A)}2a) 

{(^c(^),(^d(^0}+ = 0. (A.12b) 

Fin ally, substituting the Eq. ()A.8p - ()A.lip into Eq. 
(|A.2p and (|A.4p . the continuum model of the system is 



obtained as 

H = He^H^, (A.13a) 
^^^[A(x + L)-A(x)r}, (A.13b) 

iJe = y dx<^{x)jfrnHx), (A.13c) 

where the Hamiltonian density is 

_ r ivfa^dx + A(a;)(Ta; G {x) 

" [ G {x) ivfa^d^ + A{x + L)a^ \ ' 

(A.14) 

and we have unified the order parameters with A{x) = 
(x) and A{x -\- L) = A^ (x) , with a 4-component 

spinor ^{x) = [ (j)i{x) (j)2{x) (j^six) (l)4.{x) ]^ and modi- 
fied coupling constant 
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